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Motivated by recent experiments on intensity correlations of the waves transmitted through 
' disordered media, we demonstrate that the speckle pattern from disordered photonic crystal with 

incomplete band-gap represents a sensitive tool for determination the stop-band width. We establish 
the quantitative relation between this width and the angualar anisotropy of the intensity correlation 
function. 



0> I PACS numbers: 42.25.Dd, 42.30-d, 42.70.Qs 

Introduction. A wave propagating in a random medium undergoes multiple scattering and forms a complicated 
1— H intensity pattern, commonly referred to as a specie pattern. It is described in statistical terms, with the help of 
probability distributions and correlation functions. The intensity correlation function, {SI(r)SI{r')), where SI{r) is 
i-^ , the deviation from the averaged intensity at point r, contains a short range term [1], Ci(r,r'), which oscillates on 
' a scale of the wave length and exponentially decays beyond the mean free path, I. It also contains small long-range 
1^ , terms [2-4], which become dominant for \r — r'\ ^ I. Although the theory of speckle patterns was developed some 
time ago, it is only recently that the first experimental measurements of the spatial correlator Ci(r, r') were reported 
both for microwaves [5] and optical waves [6,7]. These experiments were carried out for isotropic disordered media; 
macroscopic isotropy was also assumed in the existing theories. In this isotropic situation comparison with the theory 
has enabled the authors [6] to infer the effective refractive index, which is the only relevant parameter of the medium 
■ in the absence of disorder. 

^ , The main message of this paper is that, in a medium with underlying spatial structure, the pattern of intensity 
Q ' correlations exhibits a vastly richer behavior as compared to the isotropic case. Moreover, additional features in 
O , Ci(r,r') carry quantitavive information about this structure. As an example, we consider a disordered incomplete- 
bandgap photonic crystal and demonstrate how the band-structure parameters can be extracted from the angular 
anisotropy of the correlator Ci . 
^ ' Superficially, it may seem that a wave with frequency ur, arriving from a distant source, after having been scattered 
, by many random inhomogeneities, will lose all information about the crystal band-structure. Our point, though, is 
' that Ci(r, r') is essentially a local object: It is determined by scattering events in the vicinity of the points r, r', and 
' it is not sensitive to the "prehistory", i.e., scattering events experienced by the wave before arrival to the point r. In 
other words, Ci(r, r') can serve as a "microscop" for observation of the local interference picture (on a scale smaller 
] than /), thus, revealing the band-structure of the inherent photonic crystal. Such a microscop is particularly well 
• suited for determination of the band-structure of realistic photonic crystals, since it is not affected by the long-range 
disorder. On the other hand, long-range disorder is a generic feature of realisitic crystals, such as silica-based synthetic 
opals. 

Disorder in synthetic opals. Silica-based synthetic opals, which are the fee self-assemled arrangement of almost 
monodispersed silica spheres, play a distinguished role in fabrication of photonic band gap materials, the potential 
^ ] of which for light manipulation was first appreciated in Refs. [8,9]. This is because the opals constitute a template 
O ■ subsequently infiltrated with high refraction index material, while the spheres are selectively removed (see e.g. [10-12]). 
Due to this distinguished role, the vast majority of experimental studies of light propagation in photonic crystals has 
been carried out on opals [13]. 

In opal photonic crystals, the contrast of the dielectric constant is weak, so that photonic bandgaps are not only 
; ^ I incomplete, but the corresponding stop-bands are relatively narrow. For this reason, the photonic band structure of 
[ this materials is strongly obscured by the disorder. As a result, any reliable determination of the stop-band either from 
disorder-broadened reflection maxima or from transmission minima, which are also srongly affected by the disorder, 
is highly ambiguous. In fact, in early studies, actual extraction of the stop-band width from the data was based on 
the "thumb rules" [14,15]. 

Later studies [13,16-18] have yielded a deeper insight into the microscopic origin of the disorder in opals. Namely, 
they indicated that one should distinguish three types of the disorder 

(i) Short-range disorder due to point defects and spread in the sphere diameters. This type of the disorder causes the 
mean free path [19] I ~ Ibfim. 

(ii) Stacking faults in [111] direction [16]. This disorder has a one- dimensional character, and thus might give rise to the 
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ID in-gap states [20] within the incomplete stop-band. In fact, experimentally observed deep and sharp transmission 
minimum [16] near the stop-band center has been accounted for these states. Subsequent interpretation of reflection 
and transmission [17] as well as diffraction data [18] also relied on the prominent role played by the stacking faults, 
(iii) macroscopic domains that arc ^ 50 — 100/im in size. Such domains inavoidably emerge in course of self-assembly 
of thick enough photonic crystals [13]. Due to their presence the stop-bands are strongly inhomogeneously broadened 
- a serious complication for determination of photonic band structure. 

Summarizing, experimental studies [13,16 18] have led to significant improvement of the understanding of the 
light propagation in strongly-disordered incomplete-bandgap photonic crystals. However, no quantitative information 
about the stop-band width could be inferred from these experiments. Below we demonstrate that this width naturally 
emerges in the angular dependence of the intensity correlator Ci . 

Correlator C\ in a disordered photonic crystal. The short-range correlator Ci(r,r'), is obtained by factorizing the 
product of four fields tp{v) in the average (/(r)/(r')). This leads [2] to the relation Ci{r,r') = |C^(r,r')|^, where 

'^^ " (7(r))V2(7(r'))V2 
is the field correlation function. It follows from the Bethe-Salpeter equation that 

C^{r,r') = ^ / dri(G(r,ri))(G*(r',ri)), (2) 



where {G{r,ri)) is the disorder-averaged Green's fimction. The physical meaning of Eq. (2) is that correlation 
between points r, r' is established via scattering on an intermediate impurity, at ri. The integral Eq. (2) can be 
reduced to 

C^ir,r') = -^Im{Gir,r')), (3) 

which has the following intuitive explanation: field correlation between a pair of points is due primarily to waves that 
are scattered in the vicinity of one point and arrive to the vicinity of the other point. The average amplitude to arrive 
from r' to r is {G{r,r')), so that C.,i,{r,r') should be a linear combination of {G{r,r')) and {G*{r,r')). Eq. (3) is 
the right linear combination since it reduces to unity for r = r'. Using Eq. (3) we present the field correlator as 

where r] = ccj/e^^^l, with e standing for the background dielectric constant, and tp*^^k{T), ui^^k are the photonic crystal 
eigenmodes and eigenfrequencies, respectively. They are characterized by the wave vector fe of the first Brillouin zone 
and the band index fi = ±1. 

Suppose that stop-band corresponds to the direction of propagation along the 2:-axis. Then, the anisotropy of 
C^(r, r') is expected within a narrow angular interval, 6 = arctan(z' — z)/(p' — p) ^ 1^^"^, where 7 = Slj/uj ^ 1 is 
the dimensionless frequency width of the stop-band. For such 9, the main contribution to the sum in Eq. (4) comes 
from the small domain of fe-space around the Bragg condition « Qb/^, where 2'jt/Qb is the period along z (see 
Fig. 1). Within this domain the dispersion law, a;^,fe, can be simplified as 



__=fc^ + _+^Q^y_+(^;,^__J , (5) 

where we have introduced the Bragg length. Lb = 4cj/Qb(5w, of the decay along z exactly at the Bragg condition. 
The corresponding eigenmodes, Vju,fe can be also simplified 

-^ij^k = Kn{z) exp - Qb/2)z + ik±p] , (6) 

where are the Bloch functions. It is seen from Eq. (5) that, for a given w^^fc = ui and the same k^, the values of 
k± are very different for upper p, = 1 and lower p = —1 bands. As a result, the contribution to C^ from the upper 
band is exponentially smaller. Then, keeping only p, = —1 terms, we can express the field-field correlation function as 



C^{r,r ) = 27r3e^/2a; / / ^-^ / ^™ ' 



c^U-i{z)U-i{z') f f / jj, J [^-^(^z ~ Qb/'^) + ipki_ cos( 
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Substituting Eq. (5) into the integral Eq. (7), and performing the integration over k± and (p we arrive at 
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where Hq is a Hankel function of a zero order. The dimensionless length TZ and "detuning" A are defined as 



TZ =\r-r'\QB, A 
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For OTZ-y^^^ :$> 1, i.e. |r — r'| ^ Lb, the asymptotic expansion of the Hankel function can be used, so that the integral 
Eq. (8) takes the form 



C^{r,r')=U-i{z)U-i{z') 
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where F{0, q) = ^q/A + ^(7^/^/2) A + ^Jl + q^ . Function F{9, q) has a special point at which the first and the 

second derivatives with respect to Q, for sonic value oi = 6c, arc zero. The position of this critical point, {Oc,qc), 
depends on the detuning A and can be parametrized by a complex number w as 
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w, ^, = 2y/^(|)'^', q. = l{4w'-9) 



1/2 



For A = we have w = 3^/^/2, 9^ = (33/4/2i/2)^i/2^ ^nd qe = 2^1"^. To avaluate the integral Eq. (10) for 9 close to 
the critical point, 9 ~ 9^, the function F{9, q) should be expanded up to the third order in (g — ^c) 



i^(^,3)=<Ao + ■^ + 55^-3(^«)^ 

\3/2 



(11) 



where 59 = 9 — 9c, Sq = {9/4:w){6/wY^^{q — qc), Fq — (4'u;^ — 9) /108, and characteristic angles 9t and 9(i are 
defined as 
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Substituting Eq. (11) into Eq. (10) and performing integration over Sq, we obtain 



C^{r,r') = 



27 
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where Ai is the Airy function. 

Discussion. It is easy to see that Eq. (13) has a form similar to the field distribution near the caustics, which is the 
envelope of rays reflected by a curved surface [21]. This similarity is not accidental. Indeed, near the caustics the 
light rays cross over from geometrically allowed directions to the "shadow". Eq. (13) also describes a crossover from 
the directions 9 ^ 9c ^ 7^^^ of free propagation to the angular domain 9 < 7^/^ where propagation is forbidden due 
to the stop-band, which plays the role of "shadow" in the fc-space. Similarly to the case of caustics, the behavior of 
near the crossover 9 = 9c is described by the following exponential decay at ^ < 



C^{9) ~ exp 
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9c - 9 
9d 
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which, in the case of intensity correlation serves as an envelope of the fast oscillations. For 9 > 9c this envelope 

,3/2 ] 

— 7r/4 . Overall, there are three angular scales in the correlator C^, namely, the 



oscillates itself as cos 
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critical angle, 9c, the period of the envelope, 9dlTZ?l^, and the period, 9tlTZ, of the oscillations. At spatial distances 
|r — r'l ~ Lb all three scales are of the same order ~ 7^/^. For large |r — r'\ we have 6c ^ 9d/TZ^^^ ^ 9t/TZ. 
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Numerical examples, illustrating separation of the three scales is given in Fig. 2. Outside the stop-band, 6 3> 7^^^ 
the correlator Ci saturates at the value ex exp (— |r — r'\/l), governed by the mean free path [1]. Obviously, our main 
result Eq. (13) applies in the vicinity of each of the equivalent Bragg directions. 

When large-scale disorder (domains) broadens significantly the reflectivity peak, R(0), the dependence Ci(9) with 
all fine details remains unchanged. Even if disorder is purely short-range, the intensity correlations are much more 
sensitive to the periodic background than the reflectivity. To illustrate this, in Fig. 3 we present the dependence 
Ci{9) together with R{6), which near the Bragg condition is given by 



ei/2 _ 1 + (ei/2 + 1) 


m - 






ei/2 + 1 + (ei/2 _ 1) 


m - 







where g{9) = Lb/^I + AiO'^/ej. The dependencies R{0), Ci{9) were calculated for the effective refractive index 
e^/^ = 1.7 and Lb = I- Similar to Ci, the width of R{0) is ~ 7^/^ However, R{9) is strongly smeared when Lb = I, 
whereas Ci{9) contains oscillations. Three oscillations are well pronounced for the distances |r — r'| equal to 21, and 
31, as shown in Fig. 3. Since these values are comparable to Lb, the argument of the Hankel function in Eq. (8) is 
~ 1. Therefore, Ci = |C^P was calculated directly from (8) without using the asymptotical form Eq. (10). Therefore, 
even when Lb ^ I, comparison of the saturation values of Ci{9) at two |r — r'| allows to determine the value of the 
mean free path. Then Lb and, correspondingly, the stop-band width, Suj = 2c/e^/'^LB, can be inferred from the 
oscillations in the following way. Knowledge of the oscillation period, Tr9t{w)/TZ together with frequency and mean 
free path, allows to deduce the parameter w, then the detunig A, defined by Eq. (9), and, finally. Lb = 21 {imA}. 
Conclusion. In the present paper we have demonstrated that the intensity correlations in a disordered photonic crystal 
with incomplete bandgap uncover the underlying band structure that is obscured by the disorder, so that conventional 
approaches applicable to the "clean" photonic crystals [22,23] fail. 

It is important to emphasize that the very phenomenon of intensity correlation is a disorder induced effect and 
that the function Ci{r,r') is calculated under the condition of diffusive propagation of the wave inside the sample. 
Note that, conceptually, with Ci playing the role of a microscope, there is an analogy between the correlation of 
the fluctuations and the experimental setup [24- 26], in which the source of light is located inside the crystal and the 
directionality of the outgoing light is studied. On the contrary, the long-range contributions [2,3] to the intensity 
correlator come from scattering trajectories that explore all directions inside the crystal, and, due to isotropization, 
do not carry information about the Bragg directions. 

Note in conclusion, that sensitivity of a number of disorder-induced effects, such as coherent backscattering [19,27] 
and directionality of the exiting diffusing light [28], to the periodic background has already been discussed in the 
literature. This sensitivity has also been employed [19,28] for the estimation of the stop-band width. In this regard, 
the intensity correlations constitute a much more accurate tool. Indeed, in the shape of the coherent backscattering 
cone, the properties of periodic structure enter only through the convolution with cooperon, which describes the large- 
scale diffusive motion, whereas the correlator Ci comes from short distances. Another advantage of using intensity 
correlations over traditional methods is that Ci is insensitive to the long-range disorder due to domains that mask 
the photonic band structure [13]. Earlier it was demonstrated [29] that, isolating a single domain and measuring 
reflection from this domain, results in a drastic narrowing of the reflection spectra. We emphasize, that correlation 
analysis of the speckle pattern allows to get rid of inhomogeneous broadening in a natural way. 

Finally, an attractive avenue for the future work comes from proposed inflltrating opals with nematic liquid crystals 
to manupulate the photonic band gap externally [30] . While in reflectivity measurements reorientation of the nematic 
director manifests itself in a simple shift of the optical Bragg reflection peak [31], the structure of the correlator Ci 
might reveal the effect of birefringence on the photonic bandstructure. 
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FIG. 1. Surface of constant frequency a;(fc^, fex) = u; is shown schematically for the Bragg frequency oj = cQs/'^t 
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FIG. 2. Correlator Ci as a function of dimonsionless angle 9/2 (QbLb) ^^'^ is shown for \r — r'\ = 5Lb and for two 
dimensionless "detunings", A, defined by Eq. (9); (a) A = 0.15i (i.e., frequency is equal to the Bragg frequency and Lb = 0.3?); 
(b) A = 0.5 + 0.15z. 




FIG. 3. (a) Correlator Ci as a function of dimensionless angle 9/2 (QbLb)"'^'^^ is shown for Lb = I and two distances 
\r — r'\ = 21 (curve 1) and \r — r'\ = 31 (curve 2); (b) Angular dependence of the reflectivity is shown for Lb = I and effective 
refractive index e^^^ = 1.7. 
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